Generation of two-mode entangled states by quantum reservoir engineering 
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A method for generating entangled cat states of two modes of a microwave cavity held is proposed. 
Entanglement results from the interaction of the held with a beam of atoms crossing the microwave 
resonator, giving rise to non-unitary dynamics of which the target entangled state is a hxed point. We 
analyse the robustness of the generated two-mode photonic "cat state" against dephasing and losses 
by means of numerical simulation. This proposal is an instance of quantum reservoir engineering of 
photonic systems. 
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I. INTRODUCTION 

Quantum reservoir engineering generally labels a strat- 
egy at the basis of protocols which make use of the non- 
unitary evolution of a system in order to generate robust 
quantum coherent states and dynamics pQ. The idea is 
in some respect challenging the naive expectation, that 
in order to obtain quantum coherent dynamics one shall 
warrant that the evolution is unitary at all stages. Due to 
the stochastic nature of the processes which generate the 
target dynamics, strategies based on quantum reservoir 
engineering are in general more robust against variations 
of the parameters than protocols solely based on uni- 
tary evolution [JH3]. A prominent example of quantum 
reservoir engineering is laser cooling, achieving prepara- 
tion of atoms and molecules at ultralow temperatures by 
means of an optical excitation followed by radiative de- 
cay [3]. The concept of quantum reservoir engineering 
and its application for quantum information processing 
has been formulated in Refs. 0|6], and further pursued 
in Refs. [7H5]. Proposals for quantum reservoir engineer- 
ing of quantum states in cavity quantum electrodynamics 
[I0HT5] and many-body systems pj [21 HS1 [17] have been 
recently discussed in the literature and first experimental 
realizations have been reported [TM2P] . Applications for 
quantum technologies are being pursued [2fJrl23|. 

In this article we propose a protocol based on quantum 
reservoir engineering for preparing a cavity in a highly 
nonclassical entangled "cat-like" state. This protocol is 
applicable to the experimental setup realized in [2H US] , 
which is pumped by a beam of atoms with random ar- 
rival times. In this setup the system dynamics intrinsi- 
cally stochastic due to the impossibility of controlling the 
arrival times of the atoms, but only their rate of injec- 
tion, and the finite detection efficiency. The protocol we 
discuss allows one to generate and stabilize an entangled 
state of two modes of a microwave resonator, by means of 
an effective environment constituted by the atoms. We 
show that when the internal state of the atoms enter- 
ing the cavity is suitably prepared and external classical 
fields couple the atomic transitions, then the asymptotic 




FIG. 1. A high-finesse microwave resonator is pumped by 
a beam of atoms with random arrival times. Two modes of 
the cavity are coupled to two atomic transitions, which are 
driven by external lasers while interacting with the fields. The 
fields undergo non-unitary dynamics, whose asymptotic state 
is an entangled state as in Eq. Q. These dynamics could be 
implemented in the experimental setup of Ref. |27| . 

state of the cavity modes takes the form 

IV-oo) = (HaHb + I - a) A \ - a) B )/N, (1) 

where \a)j denotes a coherent state of mode j — A.B 
with complex amplitude a and Af = \/2[l + exp(— 4|a| 2 )] 
is the normalization constant. 

Our proposal extends previous works of some of us, 
which are focussed on generating two-mode squeezing in 
a microwave cavity |10] and entangling two distant cav- 
ities using a beam of atoms [H]. The state of Eq. 
whose robust generation is proposed here is not simply 
entangled but possesses strongly nonclassical features, 
being a nonlocal macroscopic superposition state simi- 
lar to those discussed in Ref. [26] . The setup we consider 
is sketched in Fig. [T] and is similar to the one realized in 

Ref. main]. 

This work is structured as follows. In Sec. [ll]we sketch 
the general features of our proposal. Section [Tn| presents 
a method to engineer each of the target dynamics starting 
from the Hamiltonian of an atom of the beam, which 
interacts with the cavity for a finite time. Results from 
numerical simulations are reported and discussed in Sec. 
HVl The conclusions are drawn in Sec. [V] 
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II. TARGET MASTER EQUATION AND 
ASYMPTOTIC STATE 

Let p be the density matrix for the degrees of freedom 
of the two cavity modes and p^ = |'0oo)(V'oo| the target 
state we want to generate with |V>oo) in Eq. 0. The 
purpose of this section is to derive the master equation 

liP = Cp, (2) 

for which p^ is a fixed point, namely, 

C Poo = . (3) 

In order to determine the form of the Lindbladian C we 
first introduce the operators a and b which annihilate a 
photon of the cavity mode A and B, respectively. It is 
simple to show that p^ is a simultaneous right eigenop- 
erator at eigenvalue zero of the Liouvillians 

C jP = 7j (2C jP C}-{C}C j ,p}), 3 = 1,2 (4) 

with jj rates which are model-dependent and where the 
operators Cj read 

Ci = ^= -, C 2 = 2(ab-a 2 ). (5) 

In fact, \ipj) is eigenstate of C\ and C 2 with eigenvalue 0, 
CjlV'oo) = 0. The procedure we will follow aims at con- 
structing effective dynamics described by the Liouvillian 

C = d+C 2 (6) 

by making use of the interaction with a beam of atoms. 

Before we start, we shall remark on two important 
points. In first place, the state is not the unique 
solution of Eq. ^ when C = C\ + C 2 - Indeed, states 

|<*)a|o!)b an< ^ I — a ) A l — a ) B ' an< ^ an y superposition of 
these two states, are also eigenstates of both C\ and C 2 at 
eigenvalue zero. We denote the corresponding eigenspace 
by Hd, which is a subspace of the Hilbert space of all 
states of the two cavity modes. The most general sta- 
tionary state of C can be written as a statistical mixture, 
pss _ ^ d p d \\p d )(\j) d \ [28], where the sum spans over all 
the states \ipd) G H. d , & n d Pd are real and positive scalars 
such that J2dPd = 1- 

Nevertheless, for the evolution determined by the Lind- 
bladian of Eq. ^ the state poc is the unique asymptotic 
state provided that the initial state is the vacuum state 
for both cavity modes, po = \0a, 0b}{0a, 0b|. This can 
be shown using the parity operator defined as 

n+ = (-i) c + c + (7) 

with c± = (a ± b)/ \J2. Operator n + commutes with the 
operators C\ and C2, since 

d = c_ , C 2 = c\-c 2 _- 2a 2 . (8) 



Therefore, if the initial state can be written as statisti- 
cal mixture of eigenstates of II + with eigenvalue +1, the 
time-evolved state will also be a statistical mixture of 
eigenstates with eigenvalue +1, and so will be the steady 
state. In particular, \ipoo) is the only state of subspace 
Hd which is eigenstate of n + with eigenvalue +1, namely, 
n+l^oo) = IV'oo)) an d thus, under this condition, the 
asymptotic state will be pure and given by Poq. Here we 
will assume just this situation, i.e., that the cavity modes 
are initially prepared in the vacuum state, which is an 
even eigenvalue of operator IT + , and which represents a 
very natural initial condition. 

These considerations are so far applied to the ideal 
case in which the dynamics of the cavity modes density 
matrix are solely determined by Liouvillian C in Eq. (16b. 
In this article we will construct the dynamics in Eq. ^61) 
using a beam of atoms crossing with the resonator, as it 
is usual in microwave cavity quantum electrodynamics. 
We will then analyze the efficiency of generating state p^ 
at the asymptotics of the interaction of the cavity with 
the beam of atoms, taking also into account experimental 
limitations. 



III. ENGINEERING DISSIPATIVE PROCESSES 

Our starting point is the Hamiltonian for the coher- 
ent dynamics of an atom whose selected Rydberg transi- 
tions quasi-resonantly couple with the cavity modes. The 
atoms form a beam with statistical Poissonian distribu- 
tion in the arrival times. The mean velocity determines 
the average interaction time t during which each atom 
interacts with the cavity field, while the arrival rate r 
is such to warrant that rr <C 1, namely, the probability 
that two atoms interact simultaneously with the cavity 
is strongly suppressed. The master equation for the den- 
sity matrix \ describing the dynamics of the cavity modes 
coupled with one atom reads 

^- t X=^[H, X } + KlCx, (9) 

with H the Hamiltonian governing the coherent dynamics 
and 

JC X = 2a X a^ + 2b X tf - {a f a, X } - {tfb, X } (10) 

the superoperator describing decay of the cavity modes at 
rate n. The field density matrix is found after tracing out 
the atomic degrees of freedom, and formally reads p(t) = 
Tr at {x(i)}. In the following we will specify the form of 
Hamiltonian H and derive an effective master equation 
for the density matrix p of the cavity field interacting 
with a beam of atoms, which approximates the dynamics 
governed the Liouvillian C in Eq. Q. 

In the following we shall analyze separately each of the 
processes corresponding to the two types of Lindblad su- 
peroperators composing the sum in Eq. Q. Note that 
cavity losses are detrimental, as they do not preserve the 
parity n + of the state of the cavity. In the rest of this 
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section they will be neglected, their effect will be consid- 
ered when calculating numerically the efficiency of the 
protocol. 



A. Realization of the Lindblad superoperator L\. 

We now show how to implement the dynamics de- 
scribed by the Lindblad superoperator C±. For this pur- 
pose, we assume that the atomic transitions effectively 
coupling with the cavity modes form a A-type configu- 
ration of levels, as schematically represented in Fig. [2] 
The interaction of a single atom with the cavity modes 
is governed by the Hamiltonian 

H = hu) a a) a + hjjjjb^b + hu)2<J2.2 + ^303,3 (11) 
+h(g a a^ai t3 + gb^a 2 ,3 + H.c.) , 

where u a and ujb are the frequencies of the cavity modes, 
cj 2 (^3) is the energy of level |2) (|3)), here setting the 
energy of level |1) to zero, g a and gb are the vacuum 
Rabi frequencies characterizing the strength of the cou- 
pling of the dipolar transitions |1) — > |3) and |2) —> |3), 
respectively, with the corresponding cavity mode, and 
fj.fc = \j) (k\ is the spin- flip operator. In the following we 
assume that the transitions are resonant, i.e. u a = u>3 
and ujb — UJ3 — LU2 ■ 



field p at time t + t\ reads 




FIG. 2. Relevant atomic levels and couplings leading to the 
dynamics which realizes the Lindblad superoperator L\. The 
atom is prepared in state |— }, Eq. (13 1. 



In the reference frame rotating with the cavity modes, 
the Hamiltonian can be rewritten as 

Hl = ?n/2^( c t_ CT _ 3 + c % 3 + h.c.) , (12) 
9 



where g = \fg \ + g\ , c_ is defined in Eq. (|8| and er±.3 
|±)(3|, with 



9b\l) -9a\2) 
9 



ffa|l)+gb|2) 

9 



(13) 



while c' + is a superposition of modes a and b. This repre- 
sentation clearly shows that, if the atoms are injected in 
the state |— ) and interact with the resonator for a time 
ti such that gri <C 1, they may only absorb photons of 
the "odd" mode c_. More precisely, the condition to be 
fulfilled is gr\ \J iV_ + 1/2 <§; 1, where N_ is the mean 
number of photons in the odd mode, N- = (clc_). In 
this case, if p(t) is the state of the field at the instant in 
which an atom in state |— ) is injected, the state of the 



P (t + n) = p(t) + 



9a9b 2 

9 



2c_p{t)cl-{clc_ > p(t)} 



(14) 

This corresponds to the desired process, which drives the 
odd mode into the vacuum state. Here, we neglect correc- 
tions that are smaller by a factor of order g 2 T± (N- + 1/2) . 

Assuming that the atoms in state |— ) are injected at 
rate r\ with riTi <C 1, the probability of having two 
atoms simultaneously inside the cavity can be neglected. 
In this case the field evolution can be analysed on a 
coarsed-grained time scale At such that At ^ t\ and 
r\At <C 1. After expressing the differential quotient 
[p(t + At) — p(t)]/ At as a derivative with respect to time 
one recovers the master equation [TT] 



d 

dt f 



p(t) ~ 7l \2c^p(t)cl - {clc_, P{ t)} 



(15) 



which corresponds to the dynamics governed by super- 
operator C\ in Eq. Q. Here, 



71 = n 



9a9b 2 

9 



(16) 



We note that Eq. (15 1 is valid as long as higher or- 
der corrections are negligible. This condition provides 
an upper bound to the rate 71, i.e., 71 <C r^. However, it 
is not strictly necessary that the dynamics take place in 
this specific limit: One can indeed speed up the process 
of photon absorption from the odd mode taking longer 
interaction times between the atom and the cavity. In 
this case, the form of the master equation is different, 
but one could obtain absorption of photons from the odd 
mode. We refer the reader to Ref. [TI], where the re- 
quired time has been characterized for a similar proposal 
in the different regimes. 



B. Realization of the Lindblad superoperator Li. 

The dynamics described by the Lindblad operator £2, 
Eq. J6|) , can be realized using a level scheme as shown in 
We denote by u>j the frequency of the atomic state 
' 3), such that lu' 3 > uj' 2 > uj[ = 0. The transition is 
such that uj'3 = u a + lui,. 

A laser drives resonantly the transition |1') —> |3'), so 
that the frequency = lo' 3 = cu a + uj^. In the frame 
rotating at the frequency of the cavity modes the Hamil- 
tonian governing the coherent dynamics reads 




H 2 = ftACT 2 , 2 '+%>V 1 ,2'+.9fc&W3'+tte r i'3'+H.c.) 



(17) 



where A = uj 2 — uj a . We assume that 

9a \J ( n a)j9bV ( n a) ^ |A|, with (rij) the mean number 
of photons in the cavity mode j — A, B, and analyze 
the state of the cavity field after it has interacted with 



4 



9a> U a 
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|2') 



|1'> 



FIG. 3. Relevant atomic levels and couplings leading to the 
dynamics which approximates the Lindblad superoperator C,2- 
A classical field of amplitude Q drives resonantly the transi- 
tion |l'} — > |3'). This transition is also resonantly driven by 
two-photon processes, in which a photon of cavity mode A 
and a photon of cavity mode B are simultaneously absorbed 
or emitted. These dynamics dominate over one-photon pro- 
cesses by choosing the detuning |A| sufficiently larger than 
the coupling strengths g' a ,g' b . 



an atom which is injected in state The interaction 
time is denoted by r and is chosen such that |A|r 1 
and g' 2 (rij)T/\A\ <C 1. The density matrix for the cavity 
field at time t + r can be cast in the form [IT] 



P(t) 



1 ( 9'a9' b r 
8 A 



+i 



9a 



(At — sin At) [a^a, p] 



r,9a -2 

-2^2 sin 
2 [ A 



At\ 

— j (2apa) - {a f a,p}) 



[era, [a) a, p]], 



(18) 



where p(t) is the density matrix before the interaction 
and C 2 = 2(ab- a 2 ), Eq. ^. Here, a 2 = QA/{g' a g' b ), 
showing that the number of photons at the asymptotics 
is determined by Q. Equation ( |18[ ) has been derived in 
perturbation theory and by tracing out the degrees of 
freedom of the atom after the interaction. The first line 
of Eq. ( p~8[ ) describes two-photon processes leading to the 
target dynamics at a rate determined by the frequency 



72 



(0) _ 



1 



9a9b T 

A 



while the terms in the other lines are unwanted processes, 
which occur at comparable rates and therefore lead to 
significant deviations from the ideal behaviour. The sec- 
ond line of Eq. (18 1, in particular, corresponds to one- 
photon processes on the transition |1') — > |2'), leading to 
phase fluctuations of the cavity mode A. The third line 
describes losses of mode A due to one-photon processes, 
and the last line gives dephasing effects of cavity mode A 
associated with two-photon processes. Other detrimental 
processes, leading to dephasing and amplification of the 
field of cavity mode B, have been discarded under the as- 
sumption that the corresponding amplitude is of higher 



order. This assumption is correct as long as the ampli- 
tude SI, determining the number of photons, is chosen 
to be of the order of gf/A and fulfills the inequalities 
(|A|t)(J1t) > 1 and fir 4C 1. This is therefore a re- 
striction over the size of the cat state one can realize by 
means of this procedure. 

Let us now discuss possible strategies in order to com- 
pensate the effect of the unwanted terms in Eq. (18 1. 



We first consider the term in the second line. This term 
scales with g' 2 r / 'A and is larger than 72 . It can be 
compensated by means of a term of the same ampli- 
tude and opposite sign. This can be realized by con- 
sidering another atomic transition which is quasi reso- 
nant with the same cavity field, say, a third transition 
|laux) ~~ ^ |2 aU x) such that cavity mode A couples with 
strength g aux and detuning A aux with the dipolar tran- 
sition with |A aux | 3> g aux - If the atom is prepared in 
the superposition cos(y>)|l') + sin(y>)|l aux ) before being 
injected into the cavity, then the coherent dynamics are 
governed by Hamiltonian H' 2 = H% + /i aux , with 



h am = hA'a 2a 



%mx(a T CTi aux2aux + H.c.) , (19) 



which is reported apart for a global energy shift of 
the auxiliary levels. It is thus sufficient to select 
the parameters so that the condition cos 2 (ip)g' 2 / A + 
sin 2 (</j)<7 2 ux /A aux = is fulfilled, requiring that A and 
A aux have opposite signs. 

This operation does cancel part of the dephasing due 
to the dynamical Stark shift of cavity mode A. It does 
not compensate, however, the dephasing and dissipa- 
tion terms due to one-photon processes and scaling with 
5^ 2 A 2 sinAr and g' 2 /A 2 sin 2 (Ar/2), respectively. Nor 
does it cancel the term due to two-photon processes in the 



last line of Eq. (18 1, which scales with rate (g' 2 T / A) 2 /2. 



The remaining terms due to one-photon processes have a 
negligible effect for the choice of parameters we perform, 
since (g' 2 A 2 )/j^ ~ (9b T ) 2 an d wc choose g' b T 1 in 
order to warrant reasonably large rates (in other param- 
eter regimes, where this is not fulfilled, these terms could 
be set to zero by an appropriate selection of the velocity 
distribution of the injected atoms). 

The last term can be made smaller than 72°^ when 
{g' b /g' a ) 2 ~> 1. Nevertheless, this ratio cannot be in- 
creased arbitrarily, since the model we consider is valid 
as long as fir <C 1. This term can be identically canceled 
out when specific configurations can be realized, like the 
one shown in Fig. [4j In this configuration state |1') cou- 
ples simultaneously with the excited states |2') and |e) 
by absorption of a photon of mode A. The coherent dy- 
namics are now described by Hamiltonian H' = H 2 + h! 
with 



h! = hA'(T ee + Kg'a(a) a Ve + H.c.) 



(20) 



If the coupling strengths and detunings are such that 
g'a / A = —g a 2 /A', then not only the dynamical Stark 
shift cancels out, but interference in two-photon pro- 
cesses lead to the disappearance of the last line in Eq. 
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(18). Under this condition, the resulting master equa- 



tion is obtained in a coarse-grained time scale At assum- 
ing the atoms are injected in state |1') at rate r 2 with a 
velocity distribution leading to a normalized distribution 
p(r) over the interaction times t, with mean value t 2 and 
variance St such that At > r 2 + St. For r 2 At -C 1 the 
master equation reads 



d 

at* 



with 



and 



- 2 2C 2P Cl-{ciC 2 , Pf } 
-i/i [a f a, P] + /a (2apa f - {a f a, p}) 

72 = (r 2 /8)( 5 ^/A) 2 (r|+5r 2 ), 



(21) 



r-2 



fi 
h = r 2 



^ / drp(r) sin(Ar) 
* Jo 

I 2 ~At 

^ drp(r)sin 



(22) 
(23) 



When p(t) is a Dirac-5 function, namely, St — > 0, and 
r 2 A = 2nir with neN, then /i and f% vanish identically 
and the dynamics describes the target Liouville operator. 
Under the condition that St ^ 0, but e = ASt -C 2tt, 
then fi — 0(e 3 ) while f 2 = e 2 /4. In the other limit, in 
which p(r) is a flat distribution over [0,27r/A], then fi 
vanishes while f 2 — > 1/2. 




9a, U 'a 



FIG. 4. Level scheme leading to the master equation (21 1. 
The coupling to the additional level |e) allows one to can- 
cel out dephasing due to one-photon processes on transition 
|1')-H2'). 



C. Discussion 

In this section we have shown how to generate the tar- 
get dynamics by identifying atomic transitions and ini- 
tial states for which the desired multiphoton processes 
are driven. The level schemes we consider could be the 
effective transitions tailored by means of lasers. If the 
cavity modes to entangle have the same polarization but 
different frequencies, the levels which are coupled can 
be circular Rydberg states, while the coupling strengths 



<7j can be effective transition amplitudes, involving cav- 
ity and/or laser photons. The scheme then requires the 
ability to tune external fields so as to address resonantly 
two or more levels, together with the ability to prepare 
the internal state of the atoms entering the resonator. 
Depending on the initial atomic state, then, the dynam- 
ics can follow either the one described by superoperator 
C\ or C 2 . An important condition is that no more than a 
single atom is present inside the resonator, which sets the 
bound over the total injection rate, (r% +r 2 ) At <C 1. The 
other important condition is that the dynamics are faster 
than the decay rate of the cavity. For the experimental 
parameters we choose, this imposes a limit, among oth- 
ers, on the choice of the ratio gj/\A\, determining both 
the rate for reaching the ideal steady state as well as the 
mean number of photons per each mode, i.e., the size of 
the cat. 



IV. RESULTS 

We now evaluate the efficiency of the scheme, im- 
plementing the dynamics given by Eq. ^ with H = 
Hi+H' 2 , where Hi is given in Eq. ( 12 1 and H 2 = H 2 + h, 
with H 2 given in Eq. (17) while h depends on the addi- 



tional levels which are included in the dynamics in order 
to optimize it. The initial state of the cavity is the vac- 
uum, and the atoms are injected with rate r\ in state 
|1) (thus undergoing the coherent dynamics governed by 
Hi) and with rate r 2 in state |1), which depending on 
the considered scheme can be either (i) |1') when h = h', 
or (ii) cos(y)|l') -I- sin(<^)|l"), when h — h aux . The case 
h = is not reported, since the corresponding efficiency is 
significantly smaller than the one achievable in the other 
two cases. In order to determine the efficiency of the 
scheme we display the fidelity, namely, the overlap be- 
tween the density matrix and the target state \tpoc) 
as a function of the elapsed time. This is defined as 

F{t) - (^oo|Tr at { X (t)}^oo), 

where x(t) is the density matrix of the whole system, 
composed by cavity modes and atoms of the beam which 
have interacted with the cavity at time t, and Tr at de- 
notes the trace over all atomic degrees of freedom. 

For the purpose of identifying the best parameter 
regimes, we first analyze the dynamics neglecting the 
effect of cavity losses. Figure [5] displays the fidelity as 
a function of time when the dynamics are governed by 
Hamiltonian H = Hi + H' 2 for different realizations of 
H 2 and for different parameter choices, when the ampli- 
tude of the coherent state a = 1. Values of J- ~ 0.99 are 
reached when H' 2 = H 2 + h! is implemented. The fidelity 
then slowly decays due to higher order effects, which be- 
come relevant at longer times. The effect of two-photon 
processes involving mode A (which identically vanish for 
H 2 = H 2 + h') is visible in the two other curves, which 
correspond to the dynamics governed by H' 2 = H 2 + h aux 
when g' b — 10g' a (blue curve) and g' b = 3g' a (red curve). A 
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comparison between these two curves shows that detri- 
mental two-photon processes can be partially suppressed 
by choosing the coupling rate g' a sufficiently smaller than 



f/lr 




FIG. 5. Fidelity as a function of time (in units of the in- 
jection rate r = n = ^2) for a = 1, obtained by integrating 
numerically Eq. (J9J) after setting the cavity losses to zero, 
k = 0. The other parameters are g a T\ = g b r\ =0.1, g' b Ti = 
10 2 , g' b /A = 10~ 3 , Q,T 2 = 0.1. From top to bottom: The 
black curve refers to H'2 = H2 + h' with g' a = g' b , the other 
curves to H' 2 = H2 + h aux with g' b = 10g' a (blue) and g' b = 3g' a 
(red). 

Figure [6] displays in detail the optimal case where 
H' 2 = H2 + h'. The fidelity for the parameter choices 
g' b /A = 10~ 3 and g' b /A = 10~ 2 are reported, showing 
that a smaller ratio leads to larger fidelity in absence of 
cavity decay. The inset shows the corresponding fidelity 
when a = 0.5, which is notably larger: Reaching this 
target state starting from the vacuum, in fact, requires a 
shorter time, for which higher-order corrections are still 
irrelevant. 

The effect of cavity losses is accounted for in Fig. [7J 
where the full dynamics of master equation ([9| is simu- 
lated when H' 2 = H 2 + hi and for different choices of the 
ratio n/r. One clearly observes that the effect of cavity 
losses can be neglected over time scales of the order of 
10~ 2 /k, so that correspondingly larger rates 71 and 72 
are required. Considered the parameter choice, this is 
possible only by increasing the injection rate r. How- 
ever, this comes at the price of increasing the probability 
that more than one atom is simultaneously inside the 
resonator, thus giving rise to further sources of deviation 
from the ideal dynamics. 

These results show that degradation due to photon 
losses poses in general a problem to attain the target state 
([T ) : the rate of photon losses sets a maximum achievable 
fidelity, and also determines a time window during which 
the fidelity is close to the maximum, after which the en- 
tanglement is gradually lost. The effect of the photon 
losses is twofold: it leads to a decrease in the mean pho- 
ton number, and also breaks the symmetry preservation 
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FIG. 6. (a) Fidelity as a function of time (in units of the 
injection rate r — ri =T2) for a = 1, obtained by integrating 
numerically Eq. |9| after setting the cavity losses to zero, 
k = 0. The other parameters are f2r2 = 0.1, g a T\ = gbT\ = 
0.1, whereby the black curve is evaluated for g' b T2 = g' a T2 = 
10 2 and g'JA = 10~ 3 , while the red curve corresponds to 
9b T 2 — 9a T 2 = 10 and g' b /A = 10~ 2 (from top to bottom). 
The inset has been evaluated for the same parameters except 
for Qt2 = 0.05, leading to a = 0.5. 



in the evolution. The decrease in the mean photon num- 
ber can be compensated by increasing the strength O of 
the pumping in the implementation of the second Lind- 
blad operator, as long as the approximations made in 
Section [ill Bl are still valid. 



V. CONCLUDING REMARKS 

A strategy has been discussed which implements non- 
unitary dynamics for preparing a cavity in an entangled 
state. It is based on injecting a beam of atoms into a 
cavity, where the coherent interaction of the atoms with 
the cavity is a multiphoton process pumping in phase 
photons, so that the cavity modes approach asymptot- 
ically the entangled state of Eq. The procedure 
is robust against fluctuations of the number of atoms 
and interaction times. It is however sensitive against 
cavity losses: the protocol is efficient, in fact, as long 
as the time scale needed in order to realize the target 
state is faster than cavity decay. The effect of the pho- 
ton losses is twofold: it damps the mean photon number 
and also changes the parity of the state. It could be 
possible to partially revert the process by measuring the 
parity of the total photon number and then performing 
a feedback mechanism, similar to the one proposed in 
Refs. |30l 131] and which has been partially implemented 
in Refs. 32, 33 . Alternatively, one can find a dissipative 
way to stabilize a unique entangled target state without 
the need for feedback. This would require a process that 
can stabilize the parity of the photon number in the even 
mode. First studies have been performed showing some 
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increase in the final fidelity. We finally note that these 
ideas could also find application in other systems, such 
as circuit quantum electrodynamics setups [M] . 
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